Gluon-gluon elastic scattering amplitude in classical color field of colliding protons 
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We present a formalism for gluon-gluon elastic scattering in the presence of the classical color field 
of the protons in high energy collision. The classical field is obtained by solving the classical Yang- 
Mills equation in the covariant gauge and treated as a prescribed background for the quantum gluons 
involved in the scattering process. The interaction between the classical field and the quantum gluon 
modifies the gluon propagator, and, in turn, the gg — > gg amplitude. The modified gluon propagator 
is derived to the first non-zero order of the classical field using the Gaussian approximation in Color 
Glass Condensate and shown to satisfy the generalized Slavnov- Taylor identity. This formalism is 
the theoretical basis for our recently proposed classical color field modified minijet model where we 
show that the pp and pp cross section data from \fs = 5 GeV to 30 TeV can be satisfactorily fitted 
and the model predicts a (Ins) 2 behavior for large s, which saturates the asymptotic behaviour of 
Froissart bound. 

PACS numbers: 13.85.-t, 13.75.Cs, 13.60.Hb 



I. INTRODUCTION 

In this paper we present a formalism for gluon-gluon 
scattering in the presence of classical field due to the pro- 
tons in high energy pp collision. This formalism is the 
theoretical basis of our recently proposed classical field 
modified minijet model. We recall that the QCD-inspired 
minijet model was first introduced in the 70's At the 
time, it was noticed that the rise of the total cross section 
was very similar to the jet production cross section. In 
this context, it is natural to separate the total cross sec- 
tion into a soft component and a hard component. The 
hard component is to be calculated through the pQCD 
motivated minijet model. However, the minijet cross sec- 
tion rises too rapidly, namely, if one assume the gluon 
distribution as xg oc x~~ J with J > 1 for small x, the 
cross section a ~ s J ~ 1 lns ) that violates the Froissart 
bound. In order to restore unitarity, since then the mini- 
jet model has been incorporated in the eikonal model, 
usually referred to as eikonal minijet model, by various 
authors in an attempt to tame the rise and explain the 
data quantitatively @t3- 

We consider the modification of gluon-gluon scatter- 
ing amplitude of the minijet model due to interaction 
between the gluons involved in the gluon-gluon elastic 
process and the medium consists of small- a; gluon. The 
small- a; gluons are described by the classical effective the- 
ory calculated using the Gaussian approximation [8[ in 
Color Glass Condensate (CGC), while the scattering glu- 
ons see the small- a; gluon as a background field. In the 
theory of CGC, one first defines a scale in the longitudi- 
nal momentum fraction, x, or rapidity, then assumes the 
collection of the large- a; partons as a classical random 
source p which generates a classical field through Yang- 
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Mills equation. The classical field represents the field of 
the small- a; gluon. A physical observable is an average 
over the configuration space of p with a weight function 
W[p] [91-Hlj. When the scale x decreases to x' = x — dx, 
the source includes the partons from x = 1 down to 
x = x' . The new source at scale x' consists of a new 
layer of parton (mainly gluon) compared to the source at 
scale x. This new layer of source comes from the quan- 
tum gluon fluctuation and further induces a correction 
in the classical field. This corresponds to having a new 
weight function at the new scale x'. The ^-dependence 
of the weight function is governed by a renormalization 
group (RG) equation [ia4l5| . 

In our approach, we apply the idea of describing the 
small- x gluon field as a classical field and the physi- 
cal observable is an average over configurations of the 
source. However, both the classical field ans the source 
are prescribed inputs of our approach. Once longitudinal 
momentum fractions, (xi,X2), of the incident gluons in- 
volved in the scattering process are specified, the classical 
field and the course are fixed. The classical field provides 
an arena to the scattering process. The gluons in gg — > gg 
Feynman diagram do not only interact with each other 
but also with the classical field inside the protons. The 
longitudinal momentum fractions of the scattering glu- 
ons serves as the separation scale of the classical field 
and the source of their parent protons. The scattering 
amplitude will depend on the classical source of the pro- 
tons through the interaction with the classical field. The 
physical amplitude is then obtained by taking the aver- 
age over the random sources of both protons. Although 
an analytic solution of the RG equation is difficult to ob- 
tain, an approximate solution which preserves the Gaus- 
sian structure is suggested in Q (IIM model). The au- 
thors claimed that the approximation captures the rele- 
vant physics in both saturated and dilute regime: namely, 
the two-point correlation of the source [a x (x±, y±) = 
J (p a (x±,x~)p a (y±,y~))dx~dy~ obeys BFKL equation 
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for small (x± — y±) and demonstrates color neutrality at 
long range, (a:j_ — y±) ^> 1/Q s (x). The same Ansatz in 
IS] will be used in this paper to characterize the classi- 
cal field. We first calculate the modification to the gg 
amplitude up to the first non-zero order of the classical 
source and coupling constant. Through the resummation 
of the amplitude, the final amplitude is then analytically 
continued to a region where the source is strong. 

Without loss of generality, we consider only pp in the 
rest of the paper as we expect the result at high energy is 
the same as that for pp. We first solve the classical field of 
each proton to the leading order in coupling constant and 
the strength of the source assuming the source is static. 
Working in the covariant gauge, for a single proton, the 
equation of motion become Abelian [l6| so that the total 
field of the two colliding protons is the superposition of 
the solution of a signle proton. Solving the exact solution 
for two colliding hadrons is a much more difficult problem 
and the approximate solution to the next order has been 
found |17H2fj| . Nevertheless, a significant modification to 
the scattering amplitude can be found even in the leading 
order. 

With the inclusion of the classical effective theory, 
there are two types of gauge fields in the QCD La- 
grangian: the classical field and the quantum gluons 
which are involved in the 2-to-2 scattering. The cross 
term between the classical field and the quantum gluon in 
the Lagrangian introduces new interaction vertices in the 
Feynman diagram. For the quantum gluon, we choose to 
work in the background gauge which is consistent with 
the covariant gauge for the classical field, as shown in 
Section IIVI The classical field modifies only the gluon 
propagator. The modified propagator to the leading or- 
der of the classical field is derived and shown to satisfy 
the gauge invariance condition. We found that the gluon- 
gluon scattering amplitude is suppressed significantly in 
a strong background field; therefore, when the collinear 
gluon has a small x value, the minijet cross section is 
highly suppressed. This suppression is consistent with 
the possible breakdown of the collinear factorization for- 
mula in a high gluon density region. In the small- a; 
region, the cross section should be calculated with un- 
integrated gluon density in the k± factorization scheme, 
instead of the parton model. The modified amplitude 
provides the right amount of taming of the rise of minijet 
cross section and the total cross section can qualitatively 
describe the total cross section over the entire range of 
the available data, i.e. from <Js = 5 GcV to 30 TeV. A 
detail analysis of the phenomenological implication will 
be presented in our recent proposed classical color field 
modified minijet model (2~l| 

The focus of this paper is to derive the gluon-gluon 
scattering correction due to the presence of the classical 
field of the colliding protons within the Gaussian approx- 
imation of CGC. The remainder of the paper is organized 
as follows. In Sec. HH we derive the solution of the Yang- 
Mills equation of motion for the classical field of a pro- 
ton. We discuss how the observable is obtained using the 



IIM Gaussian approximation in Sec. Mil In Sec. IIVI 
the classical field correction to the gluon-gluon scatter- 
ing amplitude is derived and shows that there is only the 
modification of the gluon propagator. This work is con- 
cluded in Sec. [V] Appendix „ is devoted to discussions of 
gauge invariance in the present approach. 



II. CLASSICAL FIELD OF A PROTON 

In this section, we discuss the classical field given by 
a specific source distribution p a {x) of a proton travel- 
ling near the speed of light. Consider a proton moving 
along the positive light-cone (+z direction) with veloc- 
ity = (1, 0, 0, 1) as a color source with current density 
p. The color charge distribution inside the proton along 
the longitudinal direction is Lorentz contracted; there- 
fore, it is localized near x~~ = (Q- Due to time dilation, 
the source with larger x appears to be static as far as the 
dynamics of the small x gluon is concerned. Thus, the 
source is independent of light-cone time x + . The source 
can be written as 

J afi (x-,x ± ) =S»+p a (x-,x ± ). (1) 

The classical field A a ^ is governed by the Yang-Mills 
equation of motion 

D ab F b^ = jafl 

Despite the non-linearity of cq.([2|), a close solution can 
be obtained with the covariant gauge d^A^ 1 = and the 
static assumption. For a static source, d + J + = 0, it is 
consistent to look for a solution of A that satisfies 

A+= A- = 0. (3) 

The static condition also applies to A so that A is inde- 
pendent of x + , A^ = A^(x~ , x±); therefore, the partial 
derivative of A with respect to x + vanishes, 

<9+/P = d'A* = 0. (4) 

With eq. and (U]), the gauge condition together re- 
duces to 

d+A+ + d-A- - diAi = 

=> = 0. (5) 

Due to the finite size of the proton, we use the boundary 
condition of the field at infinity that A(|:ejJ — > oo) = 0. 
Thus, we have 

A l (x~,x 1 _)=Q. (6) 



1 Light-cone coordinate is used: x^q = -^=(x°±x 3 ) and x l LC = x 1 . 
The subscribe LC will be omitted in our notation. 
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Only A + = A- is non-zero which also implies A^A bfJ - — 
A a -A b+ + A a+ A b - - A ai A bl = 0. These conditions 
greatly simplify eq. ©. Eq. ([2} for p = + reduces to 
2D Possion equation, 



V 2 j_A a + (x-,x ± ) = p(x~,x±), 



(7) 



and the solution is 



Al + {x-,x ± ) 



(27T) 2 



d 2 k 



■(x-y)j 



(8) 

We assign an index 1 to p to indicate that the source is 
moving toward +z. For a source p 2 moving toward the 
opposite direction, the static assumption results A + = 
and d + A fl = 0. The only non-zero field is 



(2tt) 2 X A; 2 



In general, if there are two approaching sources, there will 
be both non-zero A + and A~ . Due to color precession 
of the sources, the sources cannot be treated as static 
anymore. The field of each source will induce a change 
in the other source. Since we are working toward the first 
order correction of a bare vacuum due to the source, these 
higher order cross inductions will be ignored. Therefore, 
the total field to the leading order approximation of the 
source 5»+pl + + 8 V ~ p%~ is 



A a„ = S »+ A a + [pi] + S^A a 2 -[ P2 ] + 0( PlP2 ) 



(10) 



III. 



GAUSSIAN AVERAGE 



As in the CGC approach, our source p is treated as a 
random variable. A physical observable O is calculated 
by first obtaining O = 0[p] in terms of p, then averaging 
over p with a weight functional W[p]. Using the IIM 
Gaussian Ansatz [8j , the observable O is given by 



(0)=N J VpO[p}W[p], 



where 



W[p] — exp < — / dydzi ± dz 



2X y (zi ± , z 2± ) 



(11) 



(12) 



N is normalization constant and X y (zi ±1 z 2± ) character- 
izes the correlation between two positions in the source. 

The source will transform if we switch gauge, p — > 
UpU^. The Gaussian weight function is gauge invariant 
and T>p is an invariant measure. Therefore, the observ- 
able is gauge invariant if the observable 0[p] is also in- 
variant when it is evaluated with a gauge transformed 
source p' = U pU\ therefore, 0[p] — 0[p']. We will use 
this criterion to check gauge invariance of the scattering 
amplitude in the classical field. The detail is presented 
in Appendix n 



The weight function in IIM model [8( is expressed in 
terms of variable y = ln(z~/z ) and z±. For our purpose, 
it is more convenient to use coordinate variables z = 
{z~, z±} as 



W[p] — exp <— dz dzi ± dz 2± 



p a (z ,z 1± )p a (z ,z 2± ) 



2X'(z ,zi ± ,z 2± ) 



where 



A'(0 = 



(13) 
(14) 



The function A(z~) has a support for < z-f < L~ ~ 
1/{xP + ) where xP + is the momentum of the collinear 
gluon involved in the scattering. The x value can be 
regarded as the scale that separates the source with 
x source > x and the field with Xfi e id < x. A correspond- 
ing correlation function in the transverse plane, p, can 
be defined as the integral of A over y, 



p T (zi ± - z 2± ) = 



-(x) 



dyX y (z 1± - z 2± ), 



(15) 



and its derivative relates to A(z ) as 
A 



dp y ^ dp{z~) _\_ = x'( z -). (16) 



dy 



The two-point correlation function in space-time coordi- 
nates becomes 

(p a (z 1 )p b (z 2 )} = S ab S(z^-z 2 )^(z-,z lx -z 2± ) ) (17) 

Since there is a one-to-one correspondence between the 
coordinate z~ and rapidity r = ln(l/x) and, in turn, x, 
we translate the longitudinal dependence into the rapid- 
ity dependence, namely, 

dp , - dp x 

(z ,zi± - z 2 _l) = t—{zii_ - z 2J _) (18) 



dz 



dz~ 



so that p x = r o L |fe = §^L- where L~~ - l/(xP+) 
is the longitudinal upper limit of the source. Simplifying 
the notation of the two-point correlation, we have 

(p a ( Zl )p b (z 2 )) = 5 ab 6(z^ - z 2 )fi x (z 1± - z 2± ), (19) 

where p x = p x /L~. If we integrate this correlation over 
the longitudinal directions z-f and z 2 , we recover the 
same correlation in eq. (3.5) of Q. 

The two-point correlation also provides that the Gaus- 
sian average of any term with odd power of p is zero but 
terms with even power of p is non-zero. In our calcu- 
lation, since the field A is linear to p (see eq.©), any 
non-zero leading order contribution must come from A\ 
or A\ because the sources from two different protons do 
not correlate; therefore, (^4i[pi]>l2[p2]) = 0. In terms of 
the Feynman diagram, the leading order correction corre- 
sponds to the diagram with two interaction vertices with 
the classical field. 
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IV. INTERACTION BETWEEN CLASSICAL 
AND QUANTUM GLUONS 

With the presence of the classical field due to the pro- 
ton, we write the total gluon field as the sum of the clas- 
sical field A and the quantum field B. A is treated as a 
prescribed field, a background, while B is the quantum 
gluon which is involved in scattering processes. The in- 
teractions between A and B are naturally emerged from 
the QCD Lagrangian. 

A. Lagrangian 

We write the full Lagrangian, ignoring the fermion 
field, as the classical QCD Lagrangian with A replaced 
by A + B 

C QCD = -\F a ^F* v + J a ^(A + B);, (20) 

where 

F a » v =d»(A + B) au - d u {A + B) B " 

+gf abc (A + B) bp (A + B) cu , (21) 

J is the classical source of the two protons. The La- 
grangian is invariant under the infinitesimal gauge trans- 
formation, 

{A + B)l (A> + B% = (A + B)« (22) 
+f abc {A + B)la c + ±d»a a 

provided that d^J^ 1 — which is satisfied by a static 
source. A detailed discussion of the gauge invariance of 
the theory will be presented in Section D 

We can organize the first terms of the Lagrangian in 
eq. (|2"0"|) in terms of the powers of B, 

Cgauge = -\f° F% = £ + A + £ 2 + £3 + U (23) 

where the subscripts represent the order of B contained 
in each term. Co is independent of B. Upon volume 
integration, it becomes a constant in the action, so it 
does not affect the observables and can be ignored. The 
first order term is 

d = -Df F b ^B« + (total derivative) 

= -J ap B a p + (total derivative), (24) 

according to the field equation of the classical field, eq. 
([2]), where D is the covariant derivative involving only A 

D* c = (d^ + gf abc Al) (25) 

and F is the field tensor of A, 



L\ cancels with J ■ B in the Lagrangian so that there is 
no single leg vertex diagram for B. £3 and £4 contain 
terms proportional to 

gB\ g 2 AB s , g 2 B\ 

The first type is the original QCD three-gluon vertex. 
The last two types are higher order terms in the coupling 
g. Our interest, instead, is in the leading order classical 
field effect which comes from the quadratic term as 

£2 = \b; [g^D p ac D cbp - D% c D» cb - 2 5 / acb F c ^] B b v . 

(27) 

B. Background gauge 

It has been shown that if one introduces a classical 
gauge field to the QCD Lagrangian, it is convenient to 
fix the gauge with background gauge [22j ■ We choose the 
gauge fixing function to be 

f a = D ab {A + B) bp = D ab B b ^, (28) 

where the last simplification is valid because A satisfies 
d^A^ — and A p A^ — 0; thus, the classical field also sat- 
isfies background gauage condition, D ab ' A bp = d^A^ + 

A a ^)+gf abc (A 1 +A 2 )l(A 1 +A 2 )^ = gf ab %A\+ A c 2 ~ + 

A\ + A\-) = Q. 

C. Quadratic term 

Including the background gauge fixing term 

= \BpZ c D» cb Bl + total derivative, (29) 

the quadratic term in B of the Lagrangian becomes 

£'2 =£2 + £gf 

J-Bl [g^D? c D chp - 2gf acb FC] B b v 

=2 B l \sT ( sabn + 9f acb 2A c p d p 

+ g 2 f aec f cdb A e p A dp ) - 2gf acb F Cfny ] B b . (30) 

The last three terms in eq. (|30[) are the interactions 
of interest. The term gg» u f acb B°iA c p d p B b is the ver- 
tex with one classical leg and two quantum legs. Since 
A = A\ + A 2 , this term actually represents two types of 
interactions, one with each proton. An interaction with 
only one power of A does not directly contribute because 
of the vanishing Gaussian average (A) oc (p) =0. There- 
fore, we will expect the gluon to interact with the same 
classical field through this interaction twice to give a first 
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non-zero correction at the g 2 A 2 . The next term is a ver- 
tex with four legs while two of them are classical. This 
term involves A ■ A. Since Ai ■ A\ — A2 ■ A2 — 0, it leaves 
with g 2 Ai ■ A<z- A\ and A2 belongs to two different Gaus- 
sian averaging procedures, {A\A 2 ) = (Ai)(A 2 ) = 0. So 
we need to have one more interaction with either classical 
field. However, this term is already at g 2 A 2 order. Any 
extra interaction will give one order higher correction. It 
will be ignored in our first order calculation. For the last 
term, we can write the field strength into the sum of field 
strength due to proton 1, Fi, proton 2, F 2 and the cross 
term F\ 2 . Explicitly, the cross term is 



gFff = g 2 f cde {A d > 1 A e 2 v 



A df * 



It also has a vanishing leading order contribution. The 
important consequence of not having cross term in the 
leading order is that the contributions of proton 1 and 2 
can be calculated independently. 

From here on, we will only consider the proton mov- 
ing to +z direction. The following results can be easily 
converted to the case with the proton moving to — z direc- 
tion. The classical field A has only "+" component and 
the derivative w.r.t. x + vanishes d+A + = d~A + = 0, 
so A^A bf * = and only the +i and i+ components, 
F al + = -F a+i = d l A a + , do not vanish. Eq. ([30]) is 
simplified to 



f 



C' 2 =-Bl [g^(6 ab n - gf abc 2A c+ d-) 



9 



9^9" 



c + 



Bl (31) 



The inverse of g^5 ab O in the first term gives the Feyn- 
man propagator. Let us denote the second term and the 
third term as 



■'inti 



gg^ r "c B a A + C d - B b 



(32) 



£-int 2 —91 



abc \ * 



i=l 



B a Mf-<r-frsr'-)diA 0+ & v (33) 



(34) 



There are two kinds of interactions: (1) the gluon changes 
color but not polarization and (2) the gluon changes 
both color and polarization. The Feynman rule of these 
classical-quantum vertices are for vertex (1): 

gg^r bc q -A+ c (p~q), (35) 

and for vertex (2): 

2 

- 9f abc J2^9 Vi ~ ffV~)(P< ~ 1i)A c + {p - q), (36) 

where the incoming and outgoing gluons are labeled by 
p, /i, a and q, v, 6, respectively. A +C (p — q) is the Fourier 
transform of the solution of A obtained in section [TTJ 
Since A +C (x~ ,x±) is independent of x + , A +C (p — q) cx 
5(p~ — q~). This delta function turns out to be very 
crucial to simplify the calculation as discussed below. 



D. Vertex and propagator corrections 

In pQCD, the high energy gluon-gluon elastic scatter- 
ing cross section is dominated by the t- and w-channcl 
gluon exchanges. For these exchange, the classical- 
quantum interaction introduces vertex and propagator 
corrections to the Feynman diagrams as illustrated in 
Fig. [Q Thanks to the delta function, 8(k~) of A + (k), 






FIG. 1. (a), (b) and (c) are the vertex corrections and (d) is 
the propagator correction due to the classical-quantum inter- 
action. A dotted line ending with a black dot represents the 
interaction with A + . 



A- 



.1+ 



* f 

Z\ 22 

FIG. 2. Feynman diagram of the propagator in classical field. 
Each vertex can be d„t 1 or Ci n t 2 - The dotted line represents 
the momentum exchange between the gluon and the classical 
field. 

(a), (b) and (c) do not contribute to the S- matrix. To 
see this, consider diagram [Tfc) . The amplitude is pro- 
portional to 



S(l--p'-)S(p--l--q-) 



1 



1 



l 2 +ie(p-l) 2 



-dV 



1 



oc- 



1+ 



-G-p)3 



-dV 



(37) 



2(1- -P-) 



p is on-sneii so ( = p'~ > 0. As we will see later, af- 
ter Gaussian averaging, the two classical lines are in fact 
connected with a momentum k + and k± going through. 
The momentum p, p' and q satisfy momentum conser- 
vation. Therefore, (p — q) 2 = p' 2 and it implies that 
q~ = q 2 /2p + < as q 2 < in the physical region. We 
also have l~ — p~ = —q~ > 0. So both poles in eq. (j37|) 
are in the lower half plane, provided that there is not any 
/ + in the numerator. By choosing the contour to close on 
the upper half plane, one obtains zero contribution. This 
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mathematical property applies also in diagram (b) and 
(c). Similar argument has been made in an earlier study 
on photon production of quark (23J. Therefore, there is 
no vertex correction due to the classical field. 

For diagram (d), the momentum integration of the loop 
in the propagator goes like J dl + / (l + — l 2 _/2l~ + ie/2l~) 
so it does not vanish, independent of the choice of the 
contour. 

To calculate the propagator correction, we write down 
the two-point Green's function to the leading order in 

A[p] as 

(G«l(x,y;p)) = (T{B«(x)Bl(y)}) 

+ (T{B;(x) 1 - J d i z 1 d A z 2 £ l ,Jz 1 )£ mt (z 2 )Bt(y)}) 

(38) 

where the notation (• • • ) means Gaussian averaging over 
the classical source defined in eq. (fTT|) . The interaction 
term linear in A + vanishes because (A + (x)) = 0. The 
first term is the Feynman propagator in bare vacuum. 
The second term represents Feynman diagrams as in Fig. 

m 

There are totally four different contributions coming 
from the binary pair of interactions Ci n t t and Ci n t 2 ■ We 

I 



will call them vertex 1 and 2, respectively, and denote 
the propagator correction involving interactions with two 
C-inti as 1-1 term, with two Ci n t 2 as 2-2 term and with 
the mixed vertices as 1-2 and 2-1 terms. 

E. {A a+ A b+ ) 

A crucial ingredient of the calculation is the random 
source averaging of A. Applying the IIM model, eq. (|19l) , 
on A a+ A b+ of eq. © gives 

(A a +( Zl )A b + (z2)) 

where p, x (kj_) is the Fourier transform of p, x (x± — yj_). 

F. 1-1 term 

Let us demonstrate the calculation of the propagator 
due to gluon interacting with Ci ntl twice at two positions. 
Let pi , p 2 and P3 be the momentum of the gluon lines on 
the left, middle and right, respectively, in Fig. [5] The 
Two-point function is given by 



Gtxt&v) = 2 x l y( T { J d i z^z 2 Bl{x)C inH {z 1 )C intl {z 2 )Bl{y)}) 
= -(T{J d 4 z 1 d 4 z 2 B^x)(~gf cmd B c p (z 1 )A m + (z 1 )d- B d P( Zl ))(~9f enf B e a (z2)A n + (z 2 )d- B^(z 2 ))B b l ,(y)}} 



There are four different ways to contract the fields. The 
first way is to contract 

B^x)B;( Zl ), B d "{z x )Bl{z 2 ), B^(z 2 )B b v (y), 

which gives 

= - 9 2 9„uf amd f dnb J d i z 1 d i z 2 {A m + (z 1 )A n + {z 2 )) 
d p\ d p 2 d 4 ps —i —i —i 

4 ToZM ("o_\4 TJ^J^J 



(2tt) 4 (2tt) 4 (2tt) 4 p\ p\ p\ 



x e 



(40) 

The Gaussian average of A m A n produces S mn so that 
jamdfdnb becomes f amd f dmb = -N c S ab . In order to 
carry out the integration of z\ and z 2l we write (A A) 
in 4D fourier transform as 

(A a + ( Zl )A b + (z 2 )) 

= 5 ab [ d 4 fc 1 d 4 fc 2 e- ifel ^ 1 e- 4fc2 - Z2 J F 1 (fc 1 ,fc 2 ), (41) 



where 



h 2 h 2 



F(h,k 2 ) 

= Mx(fcLJ 6(ki)5(kz )5 2 (k 1± + k 2± )S(k+ + k+) 
(2tt) 3 

With some algebraic manipulations, 

r ,+ (l)ab, \ -2 cab AT , -^2 di Pld A p2d A p 3 

Gii ^0, y) = -ig g^J N c (p 2 ) — 



(42) 



9 9 2 

PiPm 



3 -ipi -x+ip 3 -y 



(27T) 

F(pi -p 2 ,p 2 - pa)- 
(43) 



P2 and P3 combine to be (p 2 ) 2 because of the S(k~) = 
S(p2 — P3) in F which implies that the classical field does 
not carry momentum at the — component. Thus, 



Pi =P 2 =P 3 



(44) 
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Integrating over p2 , we have 

d*P2 2 1 . F(px -p 2 ,P2 -P3) 



(2tt) 3 



d P2 



2P2 Pi - Pii_ + ie J k l 



£,4 Hx{k\), 



(45) 



where we have expanded p\ in the denominator, inte- 
grated P2 with the delta function 5{pi— pj) an d replaced 
P2± — Pii. with k±. The integration can be evaluated 
by closing the contour on either the upper or the lower 
half plane (one can also calculate the principal value of 
the integration from — L to L as L — > 00), 



dpi 



ZpiPi -P 2 j 



_ 9(p-) - 8(-pi) 
2pi 



(46) 



The origin of the k± integration is from the Fourier trans- 
form of the transverse distribution of the classical source 
p (see eq. ©). When a gluon with momentum Q probes 
the transverse distribution of the source, it defines the 
smallest size that can be resolved. Therefore, the spatial 
integration of the source is from 1/Q to the radius of a 
proton, R p . In momentum space, Q 2 porives a UV-cutoff 
of the integration of the k 2 ^. The IR-cutoff is set to be 
1/Rp- Therefore, it becomes 



HQ 2 



Q 



l/Rl 



>(fcj) 
ki 



d 2 k. 



(47) 



where the explicit form of I(Q 2 ,x) depends on the k\ 
dependence of p, x . The the two-point function is 



^tx "* ^(2, y) 



a*N r 



4ir 
9^ 



I(Q 2 ,x)(8(q- 



b d g g r -ia-(x-y) 
(2tt) 4 q A 



0(~q-)) 

(48) 



Now we consider the other way to contract the B fields. 
As we mentioned above, eq. (|4"4"1) . that all the — compo- 
nents of the momentum are the same. This makes the re- 
sults of contracting the B fields in the different ways iden- 
tical. For example, if one contracts B^(x)B d p (zi) and 
B c p (zi)B%(z 2 ). The color index of f crnd becomes f dmc so 
it picks up a negative sign. However, the d~ is now acting 
on the gluon line withp! at z\, d~ e~ lpil > x ~ Zl ^ = +ip^[ . It 
is different by another negative sign. So there is no over 
all sign change between the two different contraction. As 
p^ = P2 , the two contraction become exactly the same. 
So the contribution due to the 1-1 term is 4 times of the 
result from any one of the contractions. The propagator 
correction in momentum space is given by 



-g,J ab ^I(Q 2 ,x)(e( q - 



0(-g"))- 



9 4 
(49) 



G. 1-2 and 2-1 terms 

As we have already worked out the 1-1 term, the details 
for the calculations of the rest of the terms are similar. 
We will point out a few keys in the calculation. For the 

1- 2 and 2-1 terms, The interaction term Ci n t 2 consists of 
a di on the classical field. In the 1-2 term, Cint 1 is at 
zi and £int 2 is a t %2- The derivative in Li n t 2 {z2) acts at 
position Z2 of (A(zi)A(z2)) and gives — ik±. While the 

2- 1 term with the same order of contraction will have the 
derivative acting on z\ resulting +ik±. Each contraction 
in the 1-2 term will cancel with that in the 2-1 term. 



H. 2-2 term 

The interaction vertex is symmetric under exchange of 
color index a ■<-> b together with Lorentz index fi ■(->• v. 
So there are totally 8 different contractions, including 
exchanging z\ and Zi , that give identical contributions. 
At the level of the Lorentz index, the only surviving term 

is E;j ~9 l3 9^9u- The 2 - 2 term is 



N r TT 



Gt 2 a ; u (x,y)^-^ r g;g-S 



-„-X"b f d4 Pl -ip x .{ x -y) 



1PlP\ 



Px{k\) _> 
k 2 d KA 



d 2 k. 



(50) 



Therefore, the propagator correction in momentum space 



g+22 tM = — s;.9^ ab ^^(%-) - <?(-?-)), 



q <r 



(51) 



where I'{Q 2 ,x) = J ^^d 2 k^. 



V. RESULT AND DISCUSSION 

We have found the leading order correction of the 
gluon-gluon t and w-channel scattering amplitude due to 
the classical color field in the background gauge only con- 
tributes to the propagator. The classical field introduces 
two type of interactions to the gluon. As for the propa- 
gator, single interaction with the classical field does not 
contribute because of the color neutrality assumption, 
namely the overall average of color field should be zero. 
However, the fluctuation can be non-zero so that second 
order terms contribute. Among all the second order dia- 
grams, only the 1-1 and 2-2 terms survive. The 1-1 term 
is diagonal in both color and Lorentz structure, while the 
2-2 term contains g^g^ ■ When one tries to sum a series 
of 1-1 and 2-2 terms, any series with more than one 2-2 
term will be zero. It is because when one connect these 
diagram with a bare propagator, g vp is inserted between 
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two diagrams. Therefore, a 2-2 term connecting with a 
2-2 term has a form of 



1% x 9 vp x 9 P 9\ = 5 M 9 9\ = since g~ 



-0. (52) 



Furthermore, under the eikonal approximation at high 



p, a, fj, 



(p-q), b, v 



Nevertheless, the modified amplitude can be applied 
to the minijet cross section and provides interesting phe- 
nomenological implications. After resumming the lead- 
ing order correction, a suppression factor is introduced 
in the gg — > gg amplitude. This suppression factor de- 
pends on the x values of the incident gluons of the minijet 
model. The amplitude receives a stronger suppression as 
x is smaller. We calculated the minijet cross section in 
21] and found that this modified minijet model provides 
a satisfactory fit to the total cross section of pp and pp 
for a wide range of energy and has a (Ins) 2 behaviour at 
high energy. 



a, c, A 

FIG. 3. Three gluon vertex. 

energy, the three-gluon vertex shown in Fig. [3] becomes 

2gg^p x r b . (53) 

Contracting the 2-2 term to the gluon line with momen- 
tum q forces p x to be p~ and the vertex vanishes because 
p~ = for an incoming gluon moving along the posi- 
tive light-cone. Therefore, when one applies the mod- 
ified propagator to high energy gluon-gluon scattering, 
one could ignore the 2-2 term and keep only the 1-1 term. 

The difference between the protons moving at the +z 
and the —z direction is that for the one moving at the 
— z direction, the A field is A^ = A" instead of A + . To 
obtain the interactions, one just needs to exchange the 
index + •<-> — in any fields and derivatives from the inter- 
action term of the +z case. We denote the propagators 
of the case for +z and — z as G + and G~ , respectively. 

The classical field modified gluon propagator to the 
first leading order to the field A = A^ + is 

Gd fM -Go f u (q) + G+ %{q) + G n fM 

+Gi 2 fM) + G^tM), (54) 
where Go ^ v {q) = -^g IJl vS ab is the bare propagator in vac- 



and, 



-9tj.v$' 



ab 



a s N r 



I(Q 2 ,x±)(9(q T )-9(-q T ))^. 

(55) 



Gf 2 2* (?) = ^flM ^i^W ) - °(-<l T )), 



q^q 4 



(56) 

where X-\- and x_ are the x value of the source moving 
in the +z and — z direction, respectively. 

So far, this is only the first non-zero order correction 
to the amplitude in the classical field. Higher order cal- 
culation requires a higher order solution of the classical 
field from the Yang-Mills equation and also the inclusion 
of the higher order terms in the Lagrangian. We defer 
the higher order consideration to future study. 



Appendix: Gauge invariance 

In this Appendix, we will use the functional method 
to derive a gauge invariant condition for the two-point 
Green's function of the propagator and use it to check 
gauge invariance. We generalized the Slavnov- Taylor 
identity 0, [H[ to the presence of an non-zero back- 
ground field. 



1. Gauge transformation with background 

We start with the generating functional 
Z[ V ;p}= J VBA(f(B)) 

xexp| 2 J d A x (£ (^ + S,p) + £ GF + ^B Q ^)| 

(A.l) 

where N is a normalization constant and the gauge fixing 
function / is chosen as 

(A.2) 



such that 



f a (B) = D ab B b » 



£ GF = — 2 {D: b Bif, 



(A.3) 



Since the classical field A is chosen to be a prescribed 
field, it does not transform. As a result, B takes all the 
burden of the gauge transformation. The infinitesimal 
gauge transformation becomes 

,»C p'a/J, — ^fiab _ jabc a c^pbii 



A a . Ala _ A a 



Hence, 

Sf a (x) 
8a c (y) 
1 



(A.4) 
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[S ab OS{x - y) 



+ gf adc (2A d x ^ + 0*M£ M + B d x ^ + &*B*J6(x - y) 
+9 2 f aeb f bdc Al^A x + B x )^S(x-y)] . (A.5) 
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The subscript x is an abbreviation of the argument of the 
function, e.g. A x = A(x). We define 

M ac (x 7 y) =D% [dn hc + gf bdc (A x + B x f»] 5(x - y) 

(A.6) 



with an inverse M 1 (x,y) satisfying 



dyD% [d»5 bc S(x - y) 

+ gf bdc (A x + B x ) d »8{x - y)] M- lcd (y, z) 
= S ad S(x - z) (A.7) 



The determinant of the gauge fixing function with respect 
to the infinitesimal gauge transformation is 



A(/) = det (U|y) = det(M) det(±; 



(A.8) 



The determinant of 1/g can be absorbed into the normal- 
ization constant N. This determinant can also be written 
as a ghost term in the Largangian: 



det(M) 



VcVc exp 



d X [CaPCa 



gf abc c a (2A b d^ + d»A b + B b d" + d^B 



g 2 f aeb f bdc c a A d (A + B)^c c 



(A.9) 



Additional interactions between the ghost and the clas- 
sical field A are also introduced. The connected rt-point 
Green's function of the quantum gluon is given by taking 
n-derivative with respect to the source r/ then divided by 

I 



Z evaluating at r\ = 



(-»)' 



S n Z 



7]=0 

(A.10) 



2. Slavnov- Taylor Identity 

As we mentioned in Section IIII1 the observable has to 
be gauge invariant. The transition amplitude depends 
explicitly on the Green's functions which itself depends 
on the gauge fixing. Therefore, one needs to find out a set 
of conditions on the Green's functions to ensure gauge in- 
variance. This is equivalent to restricting the generating 
functional to be invariant under gauge transformation. 
To do that, we transform the generating functional using 
Eq. (IA.4j) such that 

Z[ mP ] -+ Z'[ mP '} =N J VB'A(f(B')) 

exp li J d 4 x (C (A + B', p') + C GF {B') + t£B' ") j 

(A.11) 

Since £q is gauge invariant, 

C (A + B',p') =£ (A + B,p), 
and the Jacobian of the transformation on T>B is 1, 
VB = VB', 

the difference between Z and Z' is due to the change in 
gauge fixing term, the determinant A(/) and the source 
term rjB. 



Z' =N J VB(A + SA) cxp^ J d 4 x^£ Q +£ GF +6£ GF +^B a ^ + ^SB a ^j[ 

S -^- +l Jd 4 x (s£ aF + ffiSB*^ exp | i J d A x (^£ + £ GF + v «B a ^ j . 



+ N / VB A 



Gauge invariance requires that 
= 8Z = N I VB A 



ex P ^ / d 4 x(/: + /: GF + <s^ 



(A.12) 



(A.13) 



for all a a (x), where 



S£ GF = - {D ab B b ») i J dyM ac (x, y)a c (y) (A.14) 



and 



V a JB a » 



V «r dc (A + B) d x »-±(d^)5 ac a c (x) 
- ±0" «c* Q ) (A.15) 
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From here on we will denote A(/) as the determinant 
of M ab (x,y)/ g. After the transformation, M becomes 
M + SM and the determinant becomes, under gauge 
transformation (IA.4I) . 

A' = det(l/s) det(M + SM) 

= det (p a ^ [d»S bc + gf bdc (A x + B x ) d » 
+ gf dc f dmn {A x + B x ) m »al 
+ f bdc (dZa d x )]6(x-y)) 

Expanding A' to first order in SM, A' = 
det(l/ 5 )det(M + SM) « det(l/p)A(l + Tr^MM- 1 ]). 
One can identify ^ = Tr{SMM~ x ). Explicitly, M and 
SM are 



M a x c =Df il [d"S bc + gf bdc (A x + B x ) d ^} 
5M ac =D a x \ [gf bdc f dmn (A x + B x ) m »a n x 



+f bdc (d^a d x )} ; 



therefore, the trace becomes^ 



Tr[SMM~ l ] 

=Tr [D% [gf bdc f dmn (A x + B x ) m »a:M- lce (x,y) 
+f bdc (d x >a d x )M- lce (x,y)] S ae 6(x-y) 

= Jdx [gf akl f lcn (d£A* tl M~ lca (x,x)) 

+g 2jakljldcjdmn A k^ A + B )™ ^ M" 1 Ca { X , x)] oQ 

(A.16) 



By letting a x =J dyM ab (x, y)~x b {y) and requiring the to- 
tal contribution from eq. (|A.14[) . (IA.15|) and (|A.16[) to the 
deviation of Z to vanish, neglecting the total derivative 
term, we obtain the generalized Slavnov- Taylor identities 



5Z = 0=\ / dx [gf akl f lcn (d%A^M- lca (x,x)) +g 2 f akl f ldc f dmn A h xfl (A + B)Z lf "M- lca (x,x)]M' lne (x,y) 



- g D eb B b / + J dx r,lf adc {A + B)* M - - (<9^) S ac 



M- lce (x,y)}Z 



(A.17) 



The equation is an abbreviation of {...}Z — 

j T>B{. . . } e i f dxC o+CGF+vB ^ 

To the leading order, the first two terms in eq. (| A. 17|) 
are neglected because they are at least one g order higher 
than the last two. Differentiating the last two terms with 
respect to r]^(z) and setting all 77 = 0, we obtain 

= [Df y B b /B™ 

+ i[gr db (A + B) d / + d»5 cb ] A/ fc - e 1 (z,y)}z| r;=0 

(A.18) 

By taking a covariant derivative the second term 

becomes delta functions according to eq. (|A.6[) leading 
to 

{5f y Dt%B b /B^ + i5* d 5{z - y)}z\ v=0 = (A.19) 

By identifying the connected two-point Green's function 
as 

G b ; u (x,y) = ±{B b »B°/}Z\ v=0 



2 Tr[f(x, y)\ = J dxdyS(x - y)f(x, y) = f dxf(x, x) 



and taking the Gaussian average of eq. (|A.19|) , it leads 
to 

{D a x ^D bd »G%{x lV - P )) = - t 5 ab S\x - y) (A.20) 

This condition is analogous to the transverality condition 
of gluon Green's function in Slavnov's original derivation 
(eq. (22) and (23) in [2J]). One can check the propagator 
derived in Section 5 satisfies eq. (|A.20[) 

3. Gauge invariance of the modified propagator 

We will show the Green's function 

Go^t + (*tuu» (A-21) 

obtained in eq. (fl9"l) . to the leading order in g 2 and 
A 2 satisfies the generalized Slavnov- Taylor identity, eq. 
(|A.20p . We ignore the term containing the interaction 
of the second type, Ci n t 2 , since it does not contribute to 
the scattering amplitude as discussed in Section [IVJ For 
simplicity, only one proton will be considered. Since the 
contributions of the cross terms involving two colliding 
sources vanishes after the Gaussian average, the proof 
can be easily generalized to the case of two protons. 
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I Pi - Pa 



\bvy 



0000000*0000000 

Pi Pi 

FIG. 4. Schematic Feynman diagram for gluon GF with single 
interaction with A. 



In the LHS of eq. (|A.20I) , the Gaussian average is taken 
after acting the covariant derivative on the Green's func- 
tion (GF). When the term linear to A in the covariant 
derivative multiplies with the first order term in A of the 
GF, the Gaussian averaged product will have non-zero 
contribution. Therefore, besides the Gaussian averaged 
terms (eq. IA.21j) obtained in Section llVl one also needs 
to include the diagram that has only one interaction with 
the classical field in G of the LHS of eq. (|A.20[) before 
taking the average. 



4. Term linear in d n t 1 

We have calculated the terms of 0(A 2 ) in Section ITVl 
To check gauge invariance, we also need another term 
that vanish only after taking the Guassian average. This 
is the term that the gluon interacts with Ci ntl once. Di- 
agrammatically, we consider a gluon enters from the left 
with momentum p\ then interacts at z and leaves with 
momentum p 2 as shown in Fig. The GF is 

G Ia Jl(x,y)=z(TB;(x) J dzL intl {z)B b v {y)) 

= -igf cmd J dz(TBl{x) g ^Bl{z)A c+ d-B d p (z)B b v {y)). 

(A.22) 

There are two way to contract the fields. The sum of the 
two contributions are 



d 4 pid 4 p 2 fa +p 2 ) 



G^y)=9^r°]d*z (2?r)8 vlvl 

A c + {z)e- ipi{x - z) - tMz - y) . (A.23) 
5. Checking gauge invariant 



Let us rewrite the LHS of eq. (|A.20I) according to the 
order of A as 

{D?*&* v G$,(x,v)) 

+ (d^gf^A^G^Hx, y) + d^f^A^G^x, y)) 
+ (d^Gf 1 f l/ (x i y)). (A.24) 

The first term is the derivative on a bare Feynman prop- 
agator. The Gaussian average provides no effect to this 



term. Straight forward evaluation gives 



6 ab q ll .ud!td v 



g^o x o y / 



d 4 q —i 



-iq(x-y) 



(2tt) 4 q 2 



cab 



d A q 
(2tt)4 



-iq(x-y) 



i5 ab 8(x - y), 



which is the RHS of eq. (|A.20jl , The second term 

{dx9f bcd A c v {y)G Iad ll (x, y) 
= g 2 f bcd f ade (A c + (y)A e + (z)) 

ft- [ r n J 4 P ld4 P 2 ^ +P2 ) i Pl(x - z )- ip2 ( z -y) 

°* J (2n)* p\p\ 

Using eq. (|41l) to evaluate the Gaussian average, we have 



(d£gf bed A°»{y)G I %(x,v) 
--S i - (p 2 ) / d p 2 



(27T) 6 



1 



dpi 



p\ J 2p+p 2 - p 2 ± + ie 



d 2 k. 



Identifying the integrals of p\ and k± as the same as of 
eq. (|46l) and (|47[) . and by changing p 2 to q, we finally 
have 

(d£gf bed A et '(y)G I %(x ) y) 
=8« b °^(9( q -)-9(- q -))l(Q 2 ,x) 



d 1 g r -iq(x-y) 



(2ir) 4 q 2 



(A.25) 



The third term of eq. (|A.24j) term has the same con- 
tribution as the first term, so the sum of the two terms 



5 ab ^{0(q-)-e(-q-))l(Q\x) 



d 4 q q 
(2tt) 4 q 2 



-iq{x-y) 



(A.26) 



Let's recall the last term in eq. (|A.24|) as the coordinate 
space representation of eq. (l4T)l) , 

Gtit&y) =-9^5 ab ^I(Q 2 7 x)(6(q-) 9(-q-)) 



d 4 q q' 



-iq(x-y) 



(2*)V ■ (A - 2?) 

Applying the derivatives d£dy provides —iq^iq v which 
contracts with g^ gives q 2 canceling one of the q 2 in the 
denominator. Therefore, we have 

= -S ab ^I(Q 2 ,x)(9(q-)-9(-q-)) 

7T 



rf4 g g p -iq(x-l 



(2tt) 4 q 2 



(A.28) 
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which exactly cancels the sum of the two previous terms 
in (|A.26p . Therefore, the generalized Slavnov- Taylor 
identity for two point GF is satisfied, to the leading order 
in A 2 and g 2 , by the modified gluon propagator obtained 
in Section ITVl 

So far, we have ignored the GF with the interaction 
term £ mt2 ~ gB^F^B" ~ gB~diA+B l in eq. ([33]) for 
the practical reason that these corrections do not con- 
tribute to the scattering amplitude. However, the au- 
thors in [26| indicated that the typical scale of the back- 



ground field tensor F goes like gA 2 . Therefore, one can 
consider £j n t 2 as a higher order term compared to Cinti ■ 
Ignoring Ci n t 2 can be valid even at the level of order 
counting. This justifies our result. 
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